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TECHNICAL NOTE 3223 


AN ANALYSIS OF SHOCK-WAVE CANCELLATION AND REFLECTION 
FOR POROUS WALLS WHICH OBEY AN EXPONENTIAL 
MASS-FLOW PRESSURE-DIFFERENCE RELATION 


By Joseph M. Spiegel and Phillips J. Tunnell 
SUMMARY 


Two-dimensional oblique shock-wave theory is used to define condi- 
tions for cancellation and reflection of shock waves from porous walls. 
An exponential relation between mass flow normal to the walls and pres- 
sure differential through the walls is assumed. A porosity factor is 
defined which uniquely determines cancellation conditions and is inde- 
pendent of the exponent of the mass-flow pressure-difference relation 
but is dependent upon the amount of wall suction. For the reflection 
case an approximate explicit solution for the reflected wave strength 
is derived and, in general, is found to be a function of the flow 
exponent, the amount of wall suction, and the porosity factor of the 
porous medium. It is pointed out that the flow across a curved three- 
dimensional shock wave can be related to two-dimensional flow, so that 
information as to the cancellation conditions for three-dimensional 
disturbances can be obtained from the analysis. 


INTRODUCTION 


Porous walls are used in transonic test sections at low supersonic 
speeds for the purpose of canceling or attenuating flow disturbances 
that ordinarily reflect from solid boundaries. These disturbances can 
originate Prom a test model or from extrinsic sources farther upstream. 
Two past reports, which deal theoretically with the subject, neglect 
two important factors which reduce their usefulness in the study of 
real flows. First, in general, these reports assumed a linear relation 
between the flow normal to the wall and the pressure difference across 
the wall, which experimentally is not usually the case. Second, they 
neglected the complicating effects of the interaction between shock 
waves and boundary layer at the tunnel walls so that the problem could 
be more easily analyzed. This, in itself, is not too serious an omission 
because large portions of the boundary layer in transonic wind tunnels 
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can be removed by porous suction, and its removal is beneficial from the 
standpoint of boundary-layer interaction effects and probably power 
requirements. However, neither report makes allowance for the applica- 
tion of wall suction which is required to remove the boundary layer nor 
other possible effects related to wall suction. 


In reference 1 charts are presented which describe the wall porosity 
required for complete shock-wave absorption in the absence of boundary 
layer. Additional charts obtained by graphical methods enable the deter- 
mination of reflected shock-wave strengths. For the general case, these 
charts are based on an assumed linear relation between pressure differ- 
ence across the wall and flow normal to the wall, and in two other specific 
cases are based on the experimental characteristics of two porous mate- 
rials having nonlinear flow relations. 


In reference 2 a linear relation for the porous wall is also assumed. 
This report differs from reference 1 in that an approximate explicit 
solution for the strength of the reflected shock wave is presented as 
well as an equation describing shock-wave cancellation conditions. 


In the present report a generalized exponential relation between 
pressure difference and normal flow is assumed which adequately describes 
porous-wall calibrations obtained experimentally at the Ames Laboratory 
and those presented in references 1 and 3. On this basis, conditions for 
oblique shock-wave cancellation and reflection are derived. All results 
have provisions for use of an arbitrary amount of wall suction. 


SYMBOLS 


a speed of sound 


B) constants depending on porous-wall characteristics 


Po" P+ 
P,-Po 


D suction pressure-difference ratio, 


E 1l+D 

K dimensionless porosity factor 

M Mach number 

sj exponents depending on porous-wall characteristics, n = = - 1 


P-Po p-P, 


P pressure coefficient, 
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p static pressure 


dynamic pressure 
R  perfect-gas constant 
T absolute temperature 


V velocity along streamline 
y velocity perpendicular to free stream 


Vn velocity normal to porous walt 
a Mach angle, stat) 


B M? -1 

y ratio of specific heats (1.4 for air) 

5 equivalent wedge angle 

0 angle of shock wave 

v angle between free stream and porous wall 


p density 
Subscripts 


o upstream of initial shock wave 

1 downstream of initial shock wave 

2 downstream of reflected shock wave 
c cancellation conditions 

i initial suction 

r reflection conditions 


t tank enclosing test section 


8 suction 


' conditions at same pressure difference as K, 
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ANALYSIS 


The purpose of this analysis is to relate the flow conditions of 
Oblique shock waves to the flow characteristics of porous walls, based 
on the assumption that no boundary layer due to the tunnel walls is 
present. This situation will exist only if the boundary layer is removed 
by porous-wall suction, and therefore suction will be considered in the 


analysis. 


It is first necessary to define the boundary condition. It is 
assumed that pressure drop through the porous medium follows the exponen- 
tial relation 


P - p, = B(evn)" (1) 


For convenience, equation (1) can be written as 


1-M 
PVn L m 1 n 


or in dimensionless form a porosity factor, K, can be defined such that 
PVq /PoBO -p 
= saa P) (2 
(p-p,)/P, 
By use of the perfect-gas law 


pv. RT 
K = _ a o (2a) 
P-P+ Y 


which is a convenient form for computing K from experimental data. 
The constants of equations (1) and (2) are related by 


il ENS 
pas B 


This definition of a porosity factor, K, was first presented in refer- 
ence 1 in a slightly different form. It can be interpreted as the 
ability of a porous medium to pass air per unit pressure difference. 
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Shock-Wave Cancellation 


For purposes of analysis, it is convenient to begin with the case 
of no suction applied to the porous walls, although this is entirely 
hypothetical since, in the absence of suction, a boundary layer will 
exist at the wall. In actuality, as the suction is increased, the con- 
ditions shown in figure 1 will be more nearly approached. 


The porosity factor, K, required for cancellation of an oblique 
shock wave can be determined by considering the flow changes through the 
impinging shock wave. For shock cancellation without suction (v=0), con- 
ditions in region (1) of figure 1 must be maintained at the wall down- 
stream of region (1), and 8, becomes zero. From obligue-shock-wave 
relations (ref. b) 


gone —— (3) 


(+) 


P, “Po Mo” tan 8, tan 0, (5) 


Be l+ tan 6, tan 0, 


Therefore, the porosity factor, K, for cancellation becomes (eqs. (3), 
(4), and (5) substituted into the definition of K) 


— Pvn/Pošo m Pyv, [Pots 


(p-p,)/po. (P,-Po)/P, 
or 


Ko _ Sot(0,-5,) (6) 
YMo 


when P = py and V=Q (no initial suction). Since 6, depends only 
on 5, and My, Ke and (Pi-Po) /Po are uniquely determined when 5, and My 
are specified. In figure 2, Ke is plotted as a function of (p, Po) /Po 
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(eqs. (6) and (5), respectively), both for constant values of deflection 
angle and Mach mumber. This type of plot is convenient because porous- 
wall calibrations can be superposed directly onto this grid of cancella- 
tion conditions as is shown in the figure. It is apparent that if the 
porosity of a particular wall is specified by equation (2), then for any 
free-stream Mach number, M,, cancellation can occur at only two values 
of 6, at the most, and generally just one. This restriction can be 
removed by the use of wall suction. 


For the case of shock-wave cancellation with wall suction, the 
definition of K becomes 


s Pov, /Polo s p, V, sin(5,+v) P, 
8 nl 
(p, -p,)/p, (p,-p,)  Poto 


or 
pavi COS V + PoYno COB 8,0.V,/b,V, Po 


RO Na E Mi (7) 
j (Bi-Bo) + (Po-P,) Poo 


Using equation (3) and the oblique-shock-wave equation 


Vi cos 6, (8) 


Vo 7 cos(0,-5,) 


equation (7) can be written as 


Ke cos V t DKg sin 0, cos 5, /sin(0,-5,) (9) 


Kes = 
L+ D 
where 
P. "Po 
and 


PoVo sin v/ooao A Poſo 


— (10) 
(Po-P4) /Po (p,-p,)/P, 
which defines the porosity factor for suction alone. In the above egua- 
tions P,v,, represents the mass outflow per unit area through the 
porous walls for a given suction pressure differential, Po" P+ * It is 
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seen that K,g from equation (9) as well as Ke from equation (6) are 
point functions, and therefore cancellation can be achieved by all wall- 
calibration curves passing through these points, irrespective of their 
exponents of flow, n. Since (from eq. (2)), 


no SD 
Keg - (BPŁ) 
Po 


n 
Kg = C (Ft) 


and. 


equation (9) can be written as 


TH sin 0, cos 8 C 
(1 + D) ds LA (11) 
sin(9,-5,) C K! 


for D>O and Ce/C < 1. Equation (11) specifies the amount of suction 
pressure differential required to achieve shock-wave cancellation. It 
can be seen from equation (11) that for a particular wall calibration 
determined by C and n, cancellation can now be obtained over a range of 
M, and 5,, provided the wall suction, D, is varied in accordance with 
this equation. For the approximate case of 5, and Y approaching zero, 
equation (11) becomes 


(12) 


Since equation (12) cannot be solved explicitly for D, the equation has 
been plotted in Pigure 3 with n as the parameter for the range of 0 
to -0.5 (m = 1 to 2). 


Shock-Wave Reflection 


In the previous section the porosity factor required for shock-wave 
cancellation with suction was specified. It was shown that for a given 
porous wall, any shock wave for which Cae < l can be canceled when 
suction is used. For EJE > 1, impinging shock waves will be reflected 
as shock waves of reduced intensity. 1t is the purpose of this section 
to determine the reflected wave strength (8, or p»-p,) as a function of 
the porous-wall characteristics with suction. For this case the porosity 
factor for reflection (K,g) can be expressed as a function of conditions 
in region (2) of figure l. 


— — — — e—— — 
ö— — — — — — — — — — 


8 NACA TN 3223 


Pan /Poto pV, sin [ (5,5) tv] Po 
(po-P,) /P, (p, -p,) + (PoP)  Po%o 


Kyg 5 


or 
— PY COB V + Povno S08(8,-8,)0,V2/PoTo Po (13) 
(p,-p,) + (PoP) PoBo 


By use of equations (3) and (8), equation (13) can be written as 


K: - sin 0 sin 6 
> + a cos V + © p = ACES | cos(5,-5,) 
Krs. Ke P+7Po Ke sin( 0-55) sin 01781) 
Ke p.-P 
1.4 Dog > 
P,”Po 
(114) 


But Ky is the porosity factor for reflection without suction and can 
be expressed as 


K = Aa sin(6, 52) no z p V, sin(B, 53) Py-Po (15) 
Polo P_"Po P-V sin 5, P_"Po 


or, using equations (3) and (8), equation (15) becomes 


in 0. sin(8,-5 
cer ska donna me REDE MNE (16) 


sin(0,-5,)sin 5, — Po "Py 


Equations (1%) and (16) can be combined to give 


sin 05 sin(5, 57) cos V + Kg D A cos(5,-5,) 
Krs  sin075z) | sine, Z — Ke sima) _ ŽI (17 
Ke p, -P 


iš Das 
1 Po 


līt is to be noted that equations (3), (%), (5), and (8) apply to 
region (2) when the appropriate subscripts are used. 


2Q 
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By use of equation (5) and appropriate subscripts, equation (17) becomes 


sin 6, sin(5,-5,) sin 6, 
Kra  sin(0,-5,) l sin 5, °° Y KE D Sino) cos(8,-8,)| nih 
Ars _ 1 


K. 


sin 8, cos 0, sin 02 


NA o TJE A PTA FRA 
sin 8, sin(0,-5,) cos(0,-5,) 


Eguation (18) expresses the porosity factor for reflection as a function 
of the reflected wave strength, 8, (9, is a function of 5, and M,). 

It is more useful, however, to ala 8. as a function of C, n, and 
the initial wave strength. By equation (2) 


£ (tem AWZ ES + (p,-Pg) + opel] 
"Ce P, "Po Ce P,7Po 


or 


a n 
Ss a Efan + A (19) 


Noting that 


SJM (20) 


the substitution of equations (19) and (20) into equation (17) gives 


sin 0, sin(8,-8,) 


— — cos Y 
C K! sin(0.-5,) sin 6, 
Ce 5 Ke n+1 
p ~ 
(a + D+ — S par _ Sin 02 _ Sin Oy cos(5,-5,) 
P,"Po sin(0,-8,) sin(0,-5,) 


(21) 


where K' is a point on the real calibration curve at the value of 
(p - p,)/p, corresponding to Ke (see fig. 2). Equation (18) becomes 


OT 


sin 9, sin(8,-8.) 


NE Ce v 
CN A By (38) 
Ce Ke 8 D+1 in 9 sin 0 

(aa EDR ds E | „DA A A co, 5) 
sin 8, sin(9,-3,) cos(0,-5,) sin(9,-8,) sin(0, -B, ) 


If the well angle, v, is sufficiently small so that cos vel, equation (22) permits 8, to be 
determined as a function of C (or K') by graphical inversion with 8,, n, and D as parameters. 
It is to be noted that C represents the calibration curve of the real wall under consideration; 
whereas Ce represents an imaginary callbration curve passing through any point K, with the 
game exponent n as the real wall. For example, let the characteristics of the porous medium 
under consideration be represented by the real well calibration curve shown in figure 2. If the 
incident wave, Koy, has a deflection angle of 1° at Mo = 1.20, then the imaginary curve is the 
dashed line passing through point Ker parallel to the real calibration curve. Points Kar and 
Kyr are calculated as examples in the appendix. 


To determine an approximate explicit solution for the reflected wave strength, it is first 
assumed that 52, (pa-p,)/(p,-p,), and AP are small compared to unity where 6, = a, + AQ. 
Then, by use of the binomial series, neglecting squared terms of small quantities, and letting 
E = 1 + D, the pressure term of equation (21) becomes 


(z + JE ga (2 + Era EST ¢ œ .) 
Py “By E P, "P> 


so that equation (21) becomes 


K! 1 n+l P -P 
— (1 + Beh Bea a (l-8, cot 5, + B,8,)cos y + 
Ko Pi-Po 


nmi K! Bin 0, cos 8, 


Ko ~sin(@,-8,) E + 3, (tan ŠI + pa) | (23) 


EZE NL VOVN 


The pressure term can be approximated by 
Po-Pa R sin 5, cos 9, sin 6, s 8, cos 0, (24) 
P,-P, sin 8, sin(0,-8,) cos(9,-6,)  B,sin(0,-8,)sin 5, 


Inserting equation (24) into equation (23) and solving for 8a, we obtain 


sin 9 
_ t n+l _ n+l L 
B, sin(6,-8,)sin 5, | eos y A E D sin(6,-8,) 50) cos 8, 


do a 
pot tan 0, sin 5 

B, cos V sin(9,-8,)(cos 5,-B, sin 8,) + K! (n+1)E*c08 8, 1-8, ~ (sin 8,+B,cos 53) | 

Ke EP (n+1) 


(25) 


Equation (25), therefore, is an approximate solution for the reflected wave strength (8) when it 
is small compared to the initial wave. As a further approximation, 5, can be assumed to approach 
zero and cos V = 1. Then equation (25) becomes 


B, E _ (ans - 942) | 


2 q o (26) 


The pressure coefficient of the reflected wave can be written to a first approximation as 
(ref. 4, p. 134) 


or by use of equations (3) and (8) 
2 sin 20 


Dyr A rake "AME 
si B, sin 2(0,-8,) 52 (27) 


zet ML VOVN 


Equations (25) and (27) can be combined to give 


Pr = 


1 Bin 9 
sin 20, sin 5, | cos v- E pała . piti — — cos s) 
Ke sin(8,-5,) 


piti 
1 : ten 8, sin 5 
c0(6,-8,) | 8, cos v sin(0,-5,)(cos 8,-B, Sin6,) + E (n+1) Ecos AL B, 22 (sin 5,+B, cos 8,) 


E2 (n+1) 


(28) 
When 8,—>0 and cos v = 1 as for equation (26), 


appe] 


m Bo ha a = (nt) | > 


With no suction (D = 0), (Ke-K')/Ke is of the same order as Pr or Ba, so that equations (26) 
and (29) become 


_ 8, (Ke-K') /Ke  B,AK/E 


> n +2 n+2 (30) 
and 
25, AK/Kę 
(me). — (31) 


In equations (25) to (31), a positive value of 8, or Py represents a compressive reflection end a 
negative value an expansive reflection. A typical value of AK is indicated in figure 2. 


Cock NL VOVN 
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In the foregoing equations allowance has been made for the use of 
suction. Since the amount of air to be removed (Povno)by suction alone 
is a predetermined factor, the suction pressure differential Po-P;, and, 
therefore, D may be obtained by using eguation (2) in the form 


1 
nti 
— (4 Zono (32) 
Po C Poo 
DISCUSSION 


For the practical application of the equations developed in this 
paper, it is necessary first to obtain experimental data on the charac- 
teristics of the porous wali under consideration; that is, determine C 
and n of equation (2). On the basis of reference 5, C and n would be 
expected to be a function of the free-stream velocity and secondary fac- 
tors such as Reynolds number. Since To is a function of free-stream 
Mach number and total temperature, K and hence C and n are easily deter- 
mined by the remaining measurable quantities of equation (28). 


In the most general form without any simplifying assumptions, equa- 
tion (22) relates the porous-wall characteristics (K, C, and n) to the 
reflected wave strength (8,). As previously stated, 55 can be deter- 
mined, for a given wall calibration and suction, by graphical inversion 
when cos Y = 1. Equations (26) and (29) give explicit solutions for 
the reflected wave strengths in their most approximate form with suction. 
When 5, in equation (22) is made equal to zero, equation (11) is 
obtained, which relates the porous-wall characteristics and incident 
wave strength, p,-p,, to the suction for cancellation. By definition, 
the porosity factor through which a wall calibration curve must pass for 
shock cancellation with wall suction is given by equation (9), where K, 
is related to the incident shock wave by equation (6). 


In reference 3 the effect of suction on the sign and magnitude of 
reflected disturbances is discussed, and it is pointed out that suction 
tends to make a porous wall appear more closed. This effect can be 
observed in figure 2, since K for the calibration curves decreases 
with increasing pressure difference. In reference 3 it was also indi- 
cated that, with a single porous wall, a selected amount of suction can 
produce cancellation over a range of incident waves and Mach numbers 
which otherwise would result in reflections. As an example, let Kar 
in figure 2 represent a point for cancellation of an incident wave with- 
out suction. According to equation (30) or (31), an expansive reflec- 
tion will result in proportion to (KcrI-K'r1)/KcII» which is negative. 
Then equation (12) allows an approximate computation of D and hence 
(Po-P4), which is the wall suction required to cancel the reflected wave. 
In the more general case, an initial amount of suction, Po~Pti» may be 


— — — — — — — — — —— 
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used which does not achieve complete cancellation of the reflected wave, 
as indicated in the figure by Kesi». Then the additional amount of suc- 
tion required to make this point coincide with Kegri is (pyy-P+)» It 
is to be noted that the point Kes; may be above the real calibration 
curve, in which case the suction must be reduced in order to achieve 
cancellation. In summary, suction translates points on the cancellation 
grid of figure 2 to the right and up, which results in a more solid 
effect of the porous medium. Calculations for these examples are given 
in the appendix, example II. 


Since the analysis of this report was based on two-dimensional rela- 
tions, the question arises as to the possibility of applying these rela- 
tions to three-dimensional flows. In references 6 and 7 it was shown 
that for a given pressure ratio and bow-wave angle, there is both a 
wedge and a cone that can produce these conditions. This simply means 
that the flow immediately behind a conical bow shock is essentially two- 
dimensional. If this result is assumed to apply to a curved three- 
dimensional bow wave, then it is possible to represent this three- 
dimensional shock disturbance by an equivalent two-dimensional disturb- 
ance at the tangent point of the waves. This is shown schematically in 
the lower left of figure 2. In practice, this procedure amounts to 
determining the wave strength near the wall of a three-dimensional body 
either by direct measurement or the method of characteristics (when 
applicable), and then selecting a two-dimensional deflection angle to 
correspond to the three-dimensional pressure rise. A result of this 
procedure by direct measurement is shown in figure 2 for a fineness 
ratio 12 body of revolution (RM 12) which amounted to 0.34-percent 
blocked area in a 5- by >-inch transonic test section. This curve was 
obtained by measuring the pressure changes, (p,-P,)/Po> through the bow 
wave near the wall and plotting these values at the appropriate Mach 
numbers in figure 2. This is equivalent to computing Ke from equa- 
tion (6) and plotting it at the proper value of (p,-p,)/p,. This pro- 
cedure when applied to more complex bodies at various angles of attack 
can provide information as to the range of wall porosities that is 
needed for approximate interference-free testing of three-dimensional 
models. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 7, 1954 


NACA TN 3223 15 


APPENDIX 


NUMERICAL EXAMPLES 
Example I.- Calculation of Coordinates for Points 
Kaz and Kyr in Figure 2 


Given: 


Mo = 1.20 
8 


= 1.0° 


= 58.6° (from oblique shock relations, ref. >) 


CD 


4 


No suction (D < 0) 


= «40 
Wall calibration, K = 0.103 (=) 


To find Kot 


From equation (6) 


From equation (5) 


Pi-Po _ 7Mo tan 8, tan 6, 1.400 x 1.440 x 0.01745 x 1.638 


z == 0.0 
Po 1+tan 5, tan 0, 1 + 0.01745 x 1.638 = 


Using the wall calibration 


T -0 „40 
K'ņ = 0.103 (2) = 0.103(0.0560)7":*% = 0.326 


O 


To find KrT! 


From oblique shock relations, M, = 1.160. Then, for D = 0, equation (28) 
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sin 20, sin 5, [(Ke-K')/Ke] 


cos(6,-5,) | ostale, 8) (cos 5,-f, sin 5,) + E (n+1)cos 0, 


0.889 x 0.0175 [(0.378 - 0.326) /0.378] 


neem :0,00329 
0.536[0.588 x0.8IH-(1.00-0.588 x0.0175) 10.862 x1.60 x0.521] 

and 

| Mo? Pr = 0.700X1. 0.00329 = 0.0033 

O 

Therefore, 

P_-P P.-P P.-P 

279 . 1 9 , £4 = 0.0560 + 0.0033 = 0.0593 

Po Po Po 

and 


P-P+ -0.40 0.40 
Kpr = 0.103 (=) = 0.103(0.0593) = 0.319 


If suction were applied to the porous walls for the above example, the 
reflection strength, Py, would become greater. 


Example II.- Calculation of Suction Pressure Differential 
Required to Achieve Cancellation 


Given: 
Mo = 1.10 
5, = 1/29 
Kerr = 0-278 


computed similar to example 1 
(D4-Po) /Po = 0-0347 


m -0 „40 
P-P 
Wall calibration, K = 0.103 > 
O 


30 


NACA TN 3223 17 


Since 5, is small, equation (12) (fig. 3) can be used to calculate 
the suction pressure differential, (Po-P+) /Po> required to attain cancel- 
lation. It is first necessary to determine Ce/C which depends on Kjr' 
From the wall calibration 


Krr = 0.103(0.0347)79:* =0.395 


Then 


and from figure 3 
D = 0.258 
Therefore, the suction pressure differential required is 


Po”Pg _ P,”Po 


D = 0.0347x0.258 = 0.0090 
Po O 


Since 
P 
— = ———— = 0.0347 + 0.0090 = 0.0437 


KesIT = 0-103 (0.0437)7°"*° = 0.360 


If an initial suction pressure differential, which does not produce 
cancellation, is imposed on the test-section walls, the additional suc- 
tion required, (p,;-p,)/p,, can be determined as follows: 


Let (Py-P+i)/Po = 0.0050 be the initial suction. Then 


P++"P+ Po Pe PoP 


The point Kasi» corresponding to the initial suction, can be computed 
as follows: 


From the wall calibration 
Kai = 0.103(0.0050)'* = 0.858 


18 


Using equation (9) for 8,—>0 


Kesi 
l + Dy - l+ 0.144 


and 


Pa Pia 2 P,"Po A Po"Pty 


Po Po O 


= 0.0347 + 0.0050 


- Kerr + DiKsi _ 0.27840.144x0.858 
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= 0.353 


= 0.0397 
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Figure 2.- Theoretical porosity factor for cancellation of two-dimensional shock waves and ifs relation 
to a porous-wall calibration of the form K= Clap fp, ” 
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Figure 3.7 Approximate relation of suction pressure difference 
jo porosity constant ratio for shock-wave cancellation. 
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